Decay of capillary wave turbulence 
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We report on the observation of freely decaying capillary wave turbulence on the surface of a 
fluid. The capillary wave turbulence spectrum decay is found to be self-similar in time with the 
same power law exponent than the one found in the stationary regime, in agreement with weak 
turbulence predictions. The amplitude of all Fourier modes are found to decrease exponentially 
with time at the same damping rate. The longest wavelengths involved in the system are shown to 
be damped by viscous surface boundary layer. These long waves play the role of an energy source 
during the decay that sustains nonlinear interactions to keep capillary waves in a wave turbulent 
state. 

PACS numbers: 47.35.-i, 05.45.-a, 47.52. 47.27.-i 



I. INTRODUCTION 

Waves on ocean surface are the most common example 
of wave turbulence. Wave turbulence studies the statisti- 
cal properties of a set of interacting waves where energy 
is transferred by nonlinear interactions from the forcing 
scales down to small scales at which energy is dissipated. 
A statistical theory of wave turbulence was developed in 
the 60 's, the so-called weak turbulence theory which ex- 
hibits such an energy transfer in out-of-equilibrium situa- 
tions [U, . This theory has been applied to almost every 
context involving nonlinear waves: astrophysical plas- 
mas, surface or internal waves in oceanography, Rossby 
waves in the atmosphere, spin waves in magnetic ma- 
terials, Kelvin wave in supcrfluid turbulence, nonlinear 
optics, etc. This theory is based on hypothesis such as 
weakly nonlinear waves, infinite systems and scale sepa- 
rations between energy source and sink, which may limit 
its applicability to real systems. In the last decade, the 
stationary regime of wave turbulence has been studied 
in laboratory experiments (see Q for a recent review), 
and have shown that the validity domain of this theory 
can be questionable. Non-stationary wave turbulence, 
i.e. once the forcing is stopped from a stationary state 
of wave turbulence, has been much less study although 
it can provide crucial informations to understand how 
energy is redistributed between modes during transients 
(e.g. during a sudden change of wind for ocean waves), 
and how stationary regime of wave turbulence is reached. 

Numerous studies have been performed on non- 
stationary turbulent systems such as three-dimensional 
turbulence in wind tunnel in rotating tanks and 
two-dimensional turbulence in: plasma [y] , soap films 
and in shallow water @, as well as in quantum turbu- 
lence . For non-stationary wave turbulence, theoretical 
p].[lot and numerical [ll|,[l2| efforts have been performed, 
whereas experiments are scarce and only concern capil- 
lary wave turbulence [Tsj , or elastic wave turbulence on a 
thin plate For elastic waves, better agreement with 
weak turbulence theory has been obtained compared to 
the stationary case where the forcing induces anisotropy 
p^ . For the decay of capillary wave turbulence, theoret- 



ical works have predicted a self-similar solution in time 
for the energy spectrum [l^, [l^ that is compatible with 
the observations on decaying capillary wave turbulence 
on the surface of lic^uid hydrogen [l^ . 



In this paper, we report on an experimental study of 
freely decaying gravity-capillary wave turbulence on the 
surface of various fluids. The frequency spectrum of the 
capillary wave amplitude is found to be self-similar dur- 
ing most of the decay. The spectrum exponent of this 
non-stationary regime is found to be close to the pre- 
dicted exponent of the stationary regime. Moreover, the 
energy within the system and all the Fourier modes are 
found to decay exponentially with time. This is unusual 
in decaying turbulent systems where e nerg y is supposed 
to decay as a power law in time [l|, [l^, Il7j |. The longest 
waves present in the system are the most energetic ones 
and are shown to be damped by viscosity on the sur- 
face, with the inextensible film condition [l8|. However, 
these latter keep enough energy during the decay to sus- 
tain capillary waves in a wave turbulent regime without 
gravity wave turbulence. The amplitude of the capillary 
spectrum and its cut-off frequency are found to decrease 
as two different power laws of the total energy contained 
in the system. Finally, freely decaying capillary wave tur- 
bulence can be seen as a capillary wave turbulence in a 
stationary regime but with a decreasing total energy and 
a decreasing inertial range with time. 



The paper is organized as follows. In Sect. [TTl the ex- 
perimental setup and protocol are described. The exper- 
imental results are shown in Sect. IIIIl The wave height 
power spectrum is described and its properties are an- 
alyzed. The evolution of the Fourier modes and of the 
cut-off frequency of the spectrum are presented during 
the decay. The origin of wave dissipation is discussed 
in Sect. IIVI whereas Sect. IVl gives the conclusions. Fi- 
nally, additional experiments in the stationary case are 
presented in Appendix in order to compare stationary 
and non- stationary regimes. 
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II. EXPERIMENTAL SETUP AND PROTOCOL 

The experimental setup is close to the one used in [l^ . 
It consists of a circular plastic vessel, 22 cm in diameter, 
filled either with mercury or water to a height /i = 25 
mm. Water and mercury are used as the working fluid to 
study the role of the kinematic viscosity on decaying wave 
turbulence. The properties of mercury (resp. water) are: 
density, p = 13.6 x 10^ kg/m^ (resp. 10^ kg/m^), kine- 
matic viscosity v = 10~^ m^/s (resp. 10~^ m^/s), and 
surface tension 7 = 0.4 N/m (resp. 0.07 N/m). Sur- 
face waves are generated by a rectangular plunging wave 
maker (13 cm in length and 3.5 cm in height) driven 
by an electromagnetic vibration exciter driven by a ran- 
dom noise (in amplitude and frequency) band-pass fil- 
tered typically between 0.1 and 5 Hz. The amplitude 
of the surface waves at a given location is measured by 
a capacitive wire gauge plunging perpendicularly to the 
fluid at rest [l^. The eigenvalues of a circular vessel 
of radius R are given in jl8| . in terms of wave- numbers 
k. By using the dispersion relation for gravity capillary 
waves cj^ = {gk + -^k^) tanhfc/i {g being the acceleration 
of gravity, and u the angular frequency), the computed 
eigenvalue frequencies of a circular vessel, corresponding 
to a wavelength close to i? = 11 cm, are respectively 
= 3.2 Hz for anti-symmetrical modes, and = 3.9 
Hz for symmetrical ones, whatever the working fluid. 

To study freely decaying wave turbulence, a typical ex- 
periment is as follows. First, surface waves are generated 
during a long enough time (> 30 s) to reach a stationary 
wave turbulence state. The forcing is then stopped at 
t = 0, and the wave amplitude is recorded during a time 
T long enough to observe the wave damping up to a still 
state. T is chosen equal to 90 s for mercury and 40 s for 
water. The experiment then is automatically iterated TV 
times to generate statistics, and results are averaged (en- 
semble average denoted by (•)). We have chosen N large 
enough to get convergence of the statistics (TV = 180 for 
mercury and N = 300 for water). To analyze the dif- 
ferent steps of the decay of the wave amplitude signal, 
the signal is considered on short-time intervals t + 5t\ 
with 5t — 3.5 s for mercury and 2 s for water (temporal 
average denoted by "). Error bars in the following re- 
sults have been computed using various values oi 5t 
with 1 < i5t < 5 s, and considering different sets of ex- 
periments in water or mercury. Such a protocol allows 
us to study the statistical properties of freely decaying 
capillary wave turbulence with good accuracy. 
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FIG. 1. (Color online) Decay of the wave amplitude, ?7(t), as 
a function of time (solid line). Random forcing is working for 
t < 0, and is stopped at t = 0. Mercury. Dashed line corre- 
sponds to an exponential fit rj{t) = |77o|e~'^^ with r = 15.5 
s and |r;o| =3 mm. Inset: (□) temporal decay of the aver- 
aged standard deviation of wave amplitude, (cr^), and ( ) 

exponential fit (cr^(f)) = {(7,,(0))e~*/'^ with r = 15.5 s and 
{cr^(O)) — 1 mm. 



envelope of the wave amplitude, rjit) is roughly fitted by 
ij{t) = I770I exp(— t/r) (see dashed hue) where r and I770I 
are two fitting parameters corresponding to the damping 
time, and to the absolute value of 77 at t = 0, i.e. when 
the forcing is stopped. Moreover, the wave amplitude 
during its decay exhibits strong fluctuations around its 
zero mean value before reaching its stationary still state. 
Inset of Fig. [T] shows the temporal evolution of the aver- 
aged rms value of the wave amplitude {(Trjit)) = 
computed during each short-time intervals [t, t -\- 5t] and 
then averaged over N experiments. This rms value of 
the wave amplitude is found to decay exponentially in 
time with the same damping time t as above, and is well 
fitted by (cr^{t)) ~ (CT^(0))e^*/'^ (see dashed line in the 
inset). Note that similar results are found for water, r 
being the characteristic damping time for a given fluid 
(r = 5 s for water and r = 15.5 s for mercury). The 
variance ((T,j(t))^ will be considered as an estimate of the 
total energy within the system. Indeed, energy of gravity 
waves , and is greater than capillary energy. As a 
consequence, the total energy is found to decay as e"'^*/'^ 
with a t/2 damping time. 



III. NON-STATIONARY WAVE TURBULENCE 
A. Temporal decay of the wave amplitude 

Figured] shows the decay of the wave amplitude ri[t) on 
the surface of mercury as a function of time for one sin- 
gle realization once the forcing is stopped at i = 0. The 



B. Time- frequency decay of the spectrum 

To study the decay of the wave amplitude simultane- 
ously in time and in frequency, we use a spectrogram 
analysis, a well-known MATLAB signal processing tool. 
Using short-time Fourier transform, the spectrum of the 
wave amplitude at a time t* of the decay is computed 
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FIG. 2. (Color online) Averaged spectrogram {Sr,{f, t)) of the 
wave amplitude during the decay as a function of frequency 
and time. Color bar is a logarithmic scale of the spectrum 
amplitude. Pixelization in time corresponds to the short time 
interval St = 3.5 s. is an eigenvalue of the vessel (see text) 
remnant of the forcing frequencies (0.1-6 Hz). 



as Sr,{f,t*) = jf,^^\{t)ii{t + t')e-^'^-^f^dt over a short 
time interval, \t* ^t* + 5t\. Iterating in time all along the 
decay, we get the full time- frequency spectrum Sn{f,t) 
or spectrogram as a function of time and frequency. This 
also gives the temporal evolution of the spectrum of one 
Fourier mode at frequency /*: Sr,{f* ,t). The spectrum 
is then averaged over N experiments, (5^(/, t)), and is 
shown in Fig. [2] A vertical section in Fig. [2] at a time 
t* thus gives the wave amplitude spectrum as a function 
of frequency at one moment of the decay {Sn{f ,t*)) . A 
horizontal section in Fig. [5] at a frequency /* thus gives 
{Sriif * ,t)), the temporal evolution of the spectrogram of 
one Fourier mode during the decay. In Fig. [21 6t = 3.5 sis 
chosen to probe the spectrum for frequencies greater than 
2/6t w 0.6 Hz. Figure [5] shows that low frequencies are 
more energetic than high ones all along the decay. It also 
shows that high frequencies are damped faster than low 
ones. Finally, the most energetic frequency (see arrow in 
Fig. [5]) is 3.1 ± 0.1 Hz. It corresponds to a fundamental 
eigenvalue of the vessel, = 3.2 Hz that is within the 
remnant forcing frequencies. 



C. Temporal decay of the power spectrum 

Figure [3] shows the temporal evolution of the power 
spectrum (S'^(/, t*)) as a function of the frequency at 
different decay times t* . It corresponds to different ver- 
tical sections of the spectrogram in Fig. ^ top curve 
corresponding to the beginning of the decay, just after 
the forcing is stopped. Time increases from top to bot- 
tom, and corresponding curves display the dynamics of 
the decay in three main phases: i) a remnant of the forc- 
ing at the early beginning of the decay, ii) a self-similar 
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FIG. 3. (Color online) Averaged power spectrum (S,,(/, f*)) 
at different times t* of the decay. From top to bottom: t* = 
2, 5.5, 9, 16, 26.5, 37, 47.5, 58, and 68.5 s. Dashed lines 
are power laws fits ~ f^^'^ for gravity (6 < / < 25 Hz - 
top curve only) and ~ / for capillary waves (25 < / < 
100 Hz). Black arrows indicate the cut-off frequencies fc of 
the capillary wave cascade. Red (light gray) dashed arrows 
indicate vessel eigenvalues: and remnant of the forcing 
frequencies (0.1 - 6 Hz), and their harmonics 2f^ and Af^ 
(see text). Inset shows the rescaled wave height power spectra 



decay lasting most of the decay time, and iii) a purely 
dissipative phase at the end. 

i) Remnant of the forcing (t* < r/3 ~ 5 s): the power 
spectrum just after the forcing is stopped (top curve) is 
similar to the the one observed in the stationary regime 
[l9l |. Two frequency-power laws are observed, one 
corresponding to the gravity wave cascade (~ y-5.3±o.2 
for 6 < / < 20 Hz), and the other to the capillary wave 
cascade (- /-2-9±o.2 20 < / < 130 Hz). Crossover 
between the regimes occurs at fgc ^ 20 Hz, while 
dissipation occurs at higher frequency (/ > 130 Hz). 

ii) Self-similar phase (r/3 < t* < 3t ~ 45 s): all 
spectra have the same shape in the capillary regime and 
exhibit {Srf{f,t*)) ~ /~" with an exponent a that does 
not depend on time. Indeed, all frequency-power law 
fits are parallel in the capillary frequency range. Figure 
|4] shows that the exponent a is roughly independent 
of time for < i < 3r with a = 2.9 ± 0.2. This 
value is in good agreement with the theoretical value 
17/6 ~ 2.8 for the exponent of the stationary capillary 
wave turbulence [if. Note that no gravity power law is 
observed during the self-similar decay of capillary wave 
turbulence. Moreover, at high frequency, the power 
law spectrum directly reaches the noise level during the 
whole self-similar phase. The intercept between the 
power law and the noise level defines a cut-off frequency 
of the spectrum and this cut-off is found to decrease 
continuously during the self-similar decay (see black 
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FIG. 4. (Color online) Exponent a of the frequency power law 
spectrum of capillary wave turbulence as a function of time. 
One has a — 2.9 ± 0.2 during most of the decay {t < 3t — 45 
s). Dashed line is the theoretical prediction for stationary 
capillary wave turbulence 



arrows in Fig. [3]). Thus the high frequencies vanish first. 

iii) Purely dissipative phase (t* > 3t). At long times, 
no power law is observed in the capillary and gravity 
frequency ranges. The spectrum has a rounded shape 
for every frequency. 

Similar qualitative results are found when water is used 
as the working fluid. For both fluids (water and mer- 
cury), one has the self-similar decay of capillary wave 
turbulence during r/3 < t < 3t with r the character- 
istic damping time of each fluid, whereas the gravity 
cascade quickly vanishes when t > r/3. The time t/3 
can be interpreted as the time needed to "forget" the 
forcing conditions, or a memory time of the forcing. Al- 
though there is a decreasing of the inertial range during 
the self-similar decay of capillary wave turbulence (see 
Fig. [3]), all spectra are found to perfectly collapse when 
rescaled by the variance of the wave amplitude fluctua- 
tions, (S'^(/, t*))/(cr^(t*))^, as shown in the inset of Fig. 
121 The capillary power law cascade is visible up to 3t. 
This time corresponds to the "life time" of the cascade 
where more than 99% of the total energy of the system 
has been dissipated. Indeed, since ct,, e~*/'^, one has 
'''^(3''')/o'^(0) = e~^. Thus, we speculate that there is not 
enough energy within the system to sustain nonlinear in- 
teractions and the cascade phenomenology breaks down. 
Moreover, since the gravity wave cascade vanishes after 
a time t/3, we believe that the energy level necessary 
to sustain gravity wave interactions is much higher than 
that of capillary waves. Finally, note that two peaks are 
visible at low frequency in aU spectra of Fig. |3] (3.1 ±0.1 
Hz and 3.8±0.1 Hz). These frequencies correspond to 
vessel eigenvalue modes, = 3.2 Hz, and /|. = 3.9 Hz 
( [l8| and Sec. |n|), with wavelengths close to the vessel 
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FIG. 5. (Color online) Cut-off frequency as a function of 
time. First set of experiments: St = 1.95 s (□), 3.5 s (o). 
Second set: 5t = 1.95 s (o). Solid line is an exponential fit 
fc ~ e~"'^^ with a = 0.38 a fitting parameter and r = 15.5 s 
the damping time for mercury. Inset: same for water. First 
set: (□), second set (o), and third set (o). 5t — 1.5 s. Solid 
line is an exponential fit is fc ~ e~°'*^'^ with a = 0.38 a fitting 
parameter and r = 5 s the damping time for water. 



radius R. Note also that harmonics of the vessel eigen- 
value 2/^ and 4/]| are also observed close to the end of 
the decay (see Fig. [S]). 

D. Cut-off frequency of the capillary spectrum 

The cut-off frequency fc of the capillary spectrum is 
shown in Fig. [5] as a function of time for different sets 
of experiments, fc is deflned when the spectrum departs 
from its power law flt [phase i)] or when it reaches the 
noise level [phase ii) and iii)] (see black arrows in Fig. 
13]). The cut-off frequency is found to decrease roughly 
exponentially with time fc ^ e""*/"^ with a = 0.38 ±0.06 
a fitting parameter, r being the characteristic damping 
time obtained previously for mercury or water. Since 
(ct^) ~ e-*/^, one has fe ~ ((t2)0.i9±o.o3 ^^^^^^ ^^le whole 
decay for both fluids. Therefore the cut-off frequency is 
a power law of the amount of energy in the system. This 
result is compatible with the one obtained in the case of 
stationary wave turbulence (see Appendix). 

E. Temporal decay of the Fourier modes 

Figure [H] shows the temporal evolution of the ampli- 
tude of the Fourier modes {S,f{f*,t)) as a function of 
time for different /*. It corresponds to different hor- 
izontal sections of the spectrogram in Fig. [2] This 
allows us to com pare the total energy in the system 
(CT^(i))2 = (^2(^) = {Sr,{f,t))df with the amount 
of energy contained in one single Fourier mode at /* 
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FIG. 6. (Color online) Decaying of the total energy {arj{t))'^ 
(□) and of the Fourier modes {Sn{f* ,t)) as a function of time. 
Each solid line corresponds to a mode at frequency /*: both 
top /h = 3.2 Hz and = 3.8 Hz, from the third top to the 
bottom /* = 13, 23 to 93 Hz with a 10 Hz step. Dashed 
lines are exponential fits ~ e"^*'''^ for r/3 < t < 3t (slow 
damping) and ~ e~^*/'^ for t < r/3 (fast damping). Vertical 
dashed lines corresponds to r/3 and 3r. 



jf,^'USr,if,t))df = {S^{f%t))5f with Sf = 0.24 Hz 
being the Fourier transform resolution. First, one ob- 
serves that the total energy in the system decays expo- 
nentially as {arj{t)y e"^*/"^ [see □ and the dashed line 
fit in Fig. [B] . Second, excluding the beginning and the 
end of the decay, all Fourier modes are found to decay 
exponentially with the same time scale t/2 within a 10% 
error bar. This is also the damping time of the total 
energy (see parallel dashed line fits). Thus, the damp- 
ing time is found to be independent of the frequency for 
interacting waves. This result stands in stark contrast 
to the one for non interacting waves where the damping 
time of linear sinusoidal waves is a decreasing function 
of the wave frequency as a consequence of the viscous 
damping. 

One might now wonder what feeds the capillary cas- 
cade during the decay? The decay of the vessel eigenvalue 
modes = 3.2 Hz and = 3.8 Hz are also displayed 
in Fig. [5] (both top sohd curves). As for any modes, 
these fundamental modes both decay exponentially with 
the same damping time t/2 as the one of the total en- 
ergy. Moreover, one also observes that a large amount of 
energy is contained in these vessel modes. One can esti- 
mate that roughly 70% of the total energy is contained 
in frequencies between 3 and 4 Hz during all the decay. 
As a consequence, long waves close to these modes play 
the role of an energy source that feed the capillary cas- 
cade during the decay, and the relaxation of these modes 
characterize the global decay of the system. The energy 
decay of these long waves leads to a decrease of the in- 
ertial range of the capillary wave turbulence (see Figs. [3] 
andO. 



Let us now focus on the decay of Fourier modes at 
the very beginning and end of the decay. Fig. IH] shows 
that the exponential decay of Fourier modes, within the 
capillary inertial range, is faster for earlier times, t < t/3, 
than for t/3 < t < 3t, and seems also to be independent 
of the frequency. The typical time ~ t/5 of this fast 
damping could be a measurement of the nonlinear time 
scale of gravity waves and is found to be smaller than 
the damping time r/2. The fast damping is also related 
to the remnant of the forcing up to a time t/3 for which 
50% of the total energy that has been dissipated. After 
this time, the gravity power law spectrum is no longer 
observed, meaning that there is not enough energy to 
sustain the gravity wave cascade anymore. The capillary 
wave turbulence cascade is observed for longer because 
less energy is needed, and this energy can come from the 
slow relaxation of the long waves. For t > 3t, less than 1 
% of the energy remains in the system, and the capillary 
wave turbulence is no longer observed in Fig. [3] 

Note that when performing similar experiments with 
more viscous liquids (aqueous solution of glycerol or sili- 
con oil), no wave turbulence regime is observed, the dis- 
sipation time is found to depend on frequency (see Sect. 
IIV[) . and is larger than the total energy damping time 
t/2. As mentioned above, this behavior is expected when 
non interacting waves are dissipated by viscosity. As a 
consequence, decaying capillary wave turbulence is not 
observed within these fiuids because the nonlinear in- 
teractions are too weak with respect to the dissipation 
process. 



IV. ORIGIN OF WAVE DISSIPATION 



Dissipation of propagating waves in a closed basin has 
been studied theoretically and experimentally by various 
authors [l^, [l^, [2l| . The wave height decay at a given fre- 
quency has been found to be exponential whatever the 
nature of the viscous dissipation. Let us define 7" the 
theoretical damping time that depends on the frequency 
and the nature of dissipation. Viscous damping in flu- 
ids can have different origins: bottom boundary layer 
(Tb), side wall boundary layer {Tw)i and surface dissi- 
pation. Two types of viscous dissipation by the surface 
are generally considered: either due to the fluid viscos- 
ity at the surface 7^~^ ^ vk'^ EM Ell j or due to a surface 
boundary layer with an incxtensible film ^ (ywY/'^k 
[isl [20j . The latter comes from the presence of surfac- 
tants / contaminants at the interface that leads to an inex- 
tensible surface where fiuid velocity should be cancelled 
at the interface. Note that these surface dissipations are 
incompatible since they corresponds to two different kine- 
matic conditions at the interface. The decay rate for the 
wave of frequency / is defined by 5 = l/(T27r/). The 
theoretical decay rate for the various types of viscous 
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where R is the size of the circular vessel, and m = 1 
the anti-symmetrical modes and m = the symmetrical 
ones. 

To understand the origin of the wave dissipation in 
our system, we have performed similar decay experiments 
with various fluids: water, mercury, silicon oils, and aque- 
ous solutions of glycerol (denoted as x%GW with x the 
glycerol percent) to vary kinematic viscosity over two or- 
ders of magnitude. Properties of these fluids are listed in 
Table m Similar wave relaxation is observed whatever the 
working fluid: the wave height decays roughly exponen- 
tially with time, and so does its rms value (cr))- The most 
energetic mode =3.2 Hz does not depend on the fluid 
used as expected, whereas the exponential decay of this 
mode gives the damping time r for each fluid as shown in 
Table HI Figure [7] shows the experimental value of decay 
rate of the frequency =3.2 Hz as a function of the fluid 
kinematic viscosity (D-symbol) , as well as the theoretical 
decay rates computed for each type of viscous dissipation 
at the same frequency. Clearly, the experimental decay 
rate scales as v^^'^ over two decades in viscosity and not 
as as expected by usual viscous dissipation. In our ex- 
periment, viscous dissipations by surface boundary layer 
and bottom boundary layer are the most important while 
the usual viscous dissipation at the surface is clearly neg- 
ligible in the damping of the fundamental mode. Friction 
at the lateral boundary is also negligible. Bottom fric- 
tion is quite important since the wavelength correspond- 
ing to is close to the vessel radius R (and thus of 
the order of 27r/i); so waves at this frequency are affected 
by the depth. Finally, the total theoretical dissipation 
6t = Ss + Sb + Sw is found in good agreement with data 
both qualitatively and quantitatively (see dashed line in 
Fig. [7|). The fact that viscous dissipation by a surface 
boundary layer takes over the usual viscous relaxation 
was previously observed in laboratory experiments with 
water [l^, [l^, [2^ , and is also realistic considering waves 
in a natural context. Indeed, if no particular attention 
is paid (such as working in clean room, filtered fluid or 
fluid with low enough surface tension), the surface dissi- 
pation by boundary layer dominates the vk^ dissipation 
[23| . Finally, note that the infinite depth condition is 
satisfied for / > 10 Hz (i.e. A < 2 cm and kh ^ 1), and 
thus bottom friction becomes also negligible for capillary 
waves. Consequently, in our experiments, the dissipation 
source for capillary waves is only due to to viscous dissi- 



Fluid 


P i^s/m') 


u (mVs) 


7 (mN/m) 


r (s) 


Mercury 


13 600 


1.1x10"^ 


400 


15.5 ±1.5 


Water 


1 000 


10"' 


73 


5±1 


20% GW 


1 020 


2x10"'' 


70 


5±1 


30% GW 


1 050 


3x10"'^ 


70 


4 ±0.5 


50% GW 


1 120 


5x10"" 


68 


3.7 ±0.4 


Silicon oil V5 


1 000 


5x10"' 


20 


4 ±0.4 


Silicon oil VIO 


1 000 


10"^ 


20 


2.5 ±0.4 



TABLE I. Physical fluid properties: density, p, kinematic vis- 
cosity, ly, and surface tension, 7 [2^ . Corresponding experi- 
mental damping time r from the temporal decay of the vessel 
mode of frequency 3.2 Hz. 
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FIG. 7. (Color online) Decay rate, l/(r27r/J|), of the wave 
mode at — 3.2 Hz as a function of viscosity f for various 
fluids (□). Theoretical decay rate by viscous dissipation: (o) 
surface boundary layer, Ss from Eq. (±) bottom surface 
layer, Sb from Eq. ((3|, (o) wall surface layer, Sw from Eq. 
(51), and (x) viscous surface 5^ from Eq. ([T]). Dashed line is 
the total theoretical dissipation St = Ss + Sb + Sw- 



pation by a surface boundary layer, and for larger waves 
mostly by surface boundary layer since a small amount 
is also dissipated within the bottom boundary layer. 



V. CONCLUSION 

We report on laboratory experiments on the decay of 
gravity-capillary wave turbulence on the surface of var- 
ious fluids. We show that the spectrum of the capillary 
wave turbulence keeps a self-similar shape during the de- 
cay, with a frequency-power law exponent close to the one 
of the stationary regime. All Fourier modes are found to 
decay exponentially in time with the same damping rate. 
A large amount of the total energy within the system is 
contained in few longest wave modes (fundamental vessel 
eigenvalues) which also decay exponentially in time due 
to viscous dissipation within a surface boundary layer. 
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FIG. 8. (Color online) Stationary regime. Power spectrum 
Sn{f) of wave height for different mean injected powers /. 
From bottom to top: 7 = 5, 7, 10, 13, 16, 19, 26 mW. Dashed 
lines indicate power law fits ~ f~^'^- Inset shows the best 
rescaled spectrum Sn/I ' . Mercury. 



FIG. 9. (Color online) Stationary regime. Cut-off frequency 
fc of the capillary wave spectrum as a function of the mean 
injected power I. Solid line is a power law fit /c ~ / . 
Mercury. 



Since a self-similar capillary wave cascade is observed in 
time, it means that the typical time of energy flux trans- 
fer through scales by nonlinear interactions is small with 
respect to the natural damping time of the waves. Such 
a self-similar cascade during the decay of capillary wave 
turbulence has been predicted theoretically How- 
ever, the total energy is predicted to decay according 
to a power law in time which is a general result 
found in various non-stationary turbulent flows [l7j . Our 
wave turbulence system thus exhibits a unusual behav- 
ior where a turbulent cascade regime is observed with a 
total energy decaying exponentially in time as a result 
of the viscous decay of the longest wave modes. It can 
be interpreted as follows: the energy necessary to sustain 
the capillary wave turbulence is small with respect to the 
energy contained in the longest modes that feed the cap- 
illary cascade. So even if the major part of the energy is 
damped by a viscous process, there is still enough energy 
cascading through the scales to maintain a capillary wave 
turbulent state provided that the nonlinear time between 
waves is shorter than the damping time. The amplitude 
of the capillary spectrum and its cut-off frequency are 
found to decrease as two different power laws of the en- 
ergy contained in the system. Both scalings are consis- 
tent with those of the stationary regime. This means that 
freely decaying capillary wave turbulence can be seen as 
a capillary wave turbulence in a stationary regime but 
with the total energy and the inertial range decreasing 
with time. 
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Appendix: Stationary wave turbulence 

Let us compare the experimental results found in the 
non-stationary case (see above) with that of the station- 
ary regime of capillary wave turbulence. To do this, addi- 
tional experiments were performed in a stationary regime 
with decreasing injected energy. The experimental setup 
in the stationary case is the same as described in Sect. [ll] 
(see also [l^ ) except that the wavemaker is continuously 
driven by a stochastic forcing (both in amplitude and 
frequency in a range of typically 0.1 to 6 Hz). The wave 
height r](t) is then recorded during 300 s, and its power 
spectrum is computed over the whole duration of the sig- 
nal recording. The force F{t) applied by the shaker to 
the wavemaker and the velocity V(t) of the wavemaker 
are measured to access to the injected power I ^ F x V 
into the waves [l9| . 

For such a stationary regime, one finds that mean in- 
jected power I ^ afy Figure |8] shows the wave height 
spectra for different forcing amplitude 5 < / < 30 mW. 
Whatever the forcing, two power laws are observed: one 
in the gravity wave regime, S,j ~ j-5.2±o.2. ^^^^^ ^^^^ -j^ 
the capillary range of frequency, S"^ ~ y-2.8±o.2 pj.g_ 
viously reported [l^. The inertial range of the capil- 
lary spectrum also increases when / is increased. The 
power law exponent of the capillary regime is found to 
be independent of /, and the inset of Fig. [8] shows that 
Sri ^ I in agreement with previous experimental 

results [H,!!!]. Since ~ T/"^ *^ and I ~ a^, one thus 
has Sjj ~ a^f~^-^ for the stationary regime which is the 
same result as the one found experimentally for freely 
decaying capillary turbulence (see Sect. IIII C|) . 

In the stationary case, the cut-off frequency fc of the 
capillary wave turbulence spectrum is defined when the 
spectrum departs from its power law fit (see Fig. [5]). 
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Figure |9] then shows that fc increases as a power law in 
7, i.e. /, - 7 ° leading to f, ~ (a2)0-2±0 05 ^^ich 

is the same result as the one found in the freely decaying 
case (see Sect. IIIID[) . 

Finally, note that in a stationary regime, the cut-off 
frequency fc is generally estimated by balancing the non- 
linear interaction scale time tml with the dissipative time 
scale td: TNLifc) ~ Toifc) [ll, H^l- However, here, we 
have no measurement of tjvl at capillary scales. On the 
other hand, the cut-off frequency of the capillary wave 
turbulence on liquid hydrogen has been reported to scale 



as fc ~ rjp with rip the wave amplitude at the forc- 
ing frequency [131 . Assuming that r]p ~ leads to 

fc ^ <^'rl^ 1 and thus a 2/3 exponent for the variance which 
differs from the 0.2 exponent found here. Several differ- 
ences between both experiments can explain this fact. In 
our case, the stochastic forcing by a wave maker leads 
to a continuous power law spectrum, whereas the para- 
metric forcing at a single frequency of Ref. [131 leads to 
a discrete spectrum of peaks of decreasing maximal am- 
plitudes. Moreover, the dissipation nature is different 
since surface viscous boundary layer is shown to be the 
main dissipative mechanism for water waves whereas a 
dissipation is assumed for liquid hydrogen [13|. 
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